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IMPACT OF Chandra CALIBRATION UNCERTAINTIES ON GALAXY CLUSTER TEMPERATURES: 

APPLICATION TO THE HUBBLE CONSTANT 

Erik D. Reese^'^, Hajime Kawahara^''', Tetsu Kitayama^, Naomi Ota"', Shin Sasaki'^ and Yasushi Suto^'"''' 



ABSTRACT 

We perform a uniform, systematic X-ray spectroscopic analysis of a sample of 38 galaxy clusters with 
three different Chandra calibrations. The temperatures change systematically between calibrations. 
Cluster temperatures change on average by roughly ~ 6% for the smallest changes and roughly ~ 13% 
for the more extreme changes between calibrations. We explore the effects of the Chandra calibration 
' on cluster spectral properties and the implications on Sunyaev-Zel'dovich effect (SZE) and X-ray 

I determinations of the Hubble constant. The Hubble parameter changes by +10% and —13% between 

O |, the current calibration and two previous Chandra calibrations, indicating that changes in the cluster 

D ' temperature basically explain the entire change in Hq. Although this work focuses on the difference 

I in spectral properties and resultant Hubble parameters between the calibrations, it is intriguing to 

note that the newer calibrations favor a lower value of the Hubble constant, Hq ~ 60 km s~^ Mpc^^, 
typical of results from SZE/X-ray distances. Both galaxy clusters themselves and the details of the 
instruments must be known precisely to enable reliable precision cosmology with clusters, which will 
be feasible with combined efforts from ongoing observations and planned missions and observatories 
covering a wide range of wavelengths. 

Subject headings: cosmic background radiation - cosmology: observations ~ distance scale - X-rays: 
. galaxies: clusters 
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^ ■ 1. INTRODUCTION 

^ • Galaxy clusters have played a major role in determining cosmological parameters. Abundances of galaxy clus- 
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ters have placed useful constraints on the huctuation 
(e.e.. IHenrv & ArnaudI Il991t IViana & Liddld Il99e 


amplitude, cr^ 
; iBahcall et al. 


, and the matter density oi the universe, 
119971 lEke et all I199SI: IBoreani et al.l I2001t 


Reior ch & BohriuEerl 120021: iSchuecker et all 120031 IHenrv 


120041: iMantz et al.l 120081: iVikhlinin et al.l 120091: iMantz et al.l 


20101 Rozo et al. 2010). Gas fraction measurements of barvonic to total mass have also placed useful constraints on 
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flM when c ombined with D/H abundance measurements of Lya clouds and b ig bang nucleosynthesis predictions (e.g.. 
White et al. 1993: David et al. 1995: White fc Fabian 1995: Neumann fc Boh ringcr 1997; Squires et al. 1 997; EvrardI 



1997HMvers'erarill997HEttori fc Fabianlll999t iMohr et al.lll999l : iGrego etl!1l200lULaRoqiie et al.TI'200^^ Such mea- 



surements provide d one of the stronges t arguments tha t non-relativistic matter alone does not close the universe 
at that time (e.g., IViana fc Liddld 119961 : lEke et al.l [1991 iKitavama fc Sutd 119961 Il997f ). Assuming the gas fraction 
' does not evolve enables another method of using gas fractions to place constraints on cosmological parameters (e.g., 
— ■ [Sasaki 1996; Pen 1997; Ettori et al. 20031: [AUen et all [2004. 2008). The determination of galaxy cluster temperatures 
is particularly impor tant because it is widely used to in fer the gravitational mass of clusters. 
In the early 1970s. iSunvaev fc Zel'dovichl (l970l[l972D discussed the inverse Compton scattering of cosmic microwave 



background (CMB) photons off of energetic electrons of the hot cluster gas causing a sma ll (< ImK) distortion 
. in the CMB spe c trum, now known as the Sunyaev-Zel'dovi ch effect (SZE) (for reviews see, iCarlstrom eTall [200l 
■ iBirkinshawl 119991 : iRephaelil [19951 : ISunvaev fc ZeldovichI [T980[ ) . The SZE is independent of redshift making it an at- 
■ ■ ' tractive tool to explore the high redshift universe. Analysis of bot h SZE and X-ray data of galaxy clusters provides 
a method of determining a direc t distance to galax y clusters (Ca valiere et al.l [19771 : iGuiiiil 119781 : [Silk fc White! 119781 : 
iCavaliere fc Fusco-Femianolll978t [iBirkinshawl I1979D . These distances are independent of the extragalactic distance 
ladder and do not rely on standard candles or rulers. SZE/X-ray distances depend only on the properties of highly 
ionized plasmas. The promise of direct distances, in large part, provided incentive to measure the small SZE signal. 

Recently, anisotropics of the CMB, the Hubble diagram of type la supernovae, and baryon acoustic oscillations 
in the galaxy power spectrum have proved to be precise indicators of cosmology. However, galaxy cluster surveys 
probe the growth of structure, one of the few methods to do so, and interpretation of survey yields has the potential to 
preci s ely constrain the equation of state of the dark energy (e.g.. lBartlett fc Silklll994lHolder et al.|[2000tlHaiman et al] 
120011 : iMajumdar fc Mohi1l2004[ ). Current and future surveys of galaxy clusters in radio (ACT, SPT, APEX-SZ; e.g.. 
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Fowler et all 120101: iHincks et al.ll2009t iFowler et alj 120071: iCarlstrom et all [200llLueker eFall 120101: iReichardt eFall 



20091) X-ray (eROSITA: e.g.. IPredehl et all 120071 [2006D . and weak-lensing (Pan-STARRS, HSC, LSST, etc.; e.g. 
Kaiseii[200l iMivazaki et al.l 120061: ll\ 



vezic et al. promise to provide unprecedented statistical samples of galaxy 

clusters that will enable precision cosmology with galaxy clusters. Therefore it is important to consider possible 

systematics inherent to using clusters as tools of cosmology. 

Traditionally SZE/X -ray estimates of the Hubbl e parameter tend t o be low (for details see, e.g.. lBirkinshawlll999l : 
ICarlstrom et al.|[2002[ ). with a few exceptions (e.g.. lMason et al.]|2001|) . However, a recent study of 38 galaxy clusters 
using data from Chandra and the OVRO/BIMA SZE imaging project find Hq 74 km s~^ Mpc~^ fBonamente et ajj 
l2006f ). consistent with o ther probes, such a s the Hubble Space Telescope (HST) Hq key project (Frccdnian ct al. 200'lX 
super nova results fe.g.. iHicken et aLl '2009'). and recent CMB primary anisotropy results (e.g., Kom atsu et al.l l2009l 
|2010[) . Previous work using the same SZE data but ROSAT and ASCA X-ray data for 18 galaxy clusters find Hq ~ 60 
km s-i Mpc-i (|Reese et al.ll2002L hereafter R02). 

We explore the effects of the Chandra calibration on cluster X-ray spectral properties and their subsequent effects on 
the inferred Hq, through an X-ray spectral analysis of Chandra archival data of SZE clusters. We perform a uniform, 
systematic spectroscopic X-ray analysis of 38 galaxy clusters, comparing the effects of three different Chandra calibra- 
tions on galaxy cluster p roperties. The clusters are chosen from a sample used for SZE/X-ray distance measurements 
(|Bonamente et all 120061 hereafter B06), facilitating the exploration of how the effects of the calibration affect the 
inferred Hubble parameter. 

B06 find much less model dependence than one might naively expect for determinati ons of ffp from SZE/X-ray 
determined distances and therefore we focus on the simplest case, the isothermal /3-model (jCavaliere fc Fusco -Femianol 
[l976, 1978). In addition, a simple isothermal model will make it easier to isolate the effects of the different Chandra 
calibrations on individual spectral properties, such as cluster gas temperature, Tg, and the X-ray cooling function. A, 
and how these differences propagate into cosmological parameter determinations. Because Hq cx T^/A{Te), the Hubble 
constant provides a good test case for tracking the effects of cluster temperature on cosmological parameters. 

The rest of the paper is organized as follows. We begin with a short pedagogical discussion on modeling of the X-ray 
emission from galaxy clusters in Section [2] with particular attention to its potential systematic uncertainties. Chandra 
data reduction and spectroscopic analysis for the three Chandra calibrations are discussed in Section [S] Implications 
for the Hubble constant are explored in Section |4] and more general implications are discussed in Section [S] Chandra 
calibrations 3.1.0, 4.1.4, and 4.2.2 are considered and are referred to more simply as 3.1, 4.1, and 4.2. Throughout 
this paper, all uncertainties are at 68% confidence and we adopt a flat, A-dominated cosmol ogy with = 0.27 
and f^A = 0.73 consistent with recent Wilkinson Microwave Anisotropy Probe (WMAP) results ()Komatsu et al.|[2010l 
|2009| ). As is the convention the Hubble constant is sometimes expressed as /i, defined as Hq = 100 h km s~^ Mpc~^. 

2. MODELING X-RAY SURFACE BRIGHTNESS OF GALAXY CLUSTERS FOR ESTIMATING THE HUBBLE 

CONSTANT 

2.1. Theoretical modeling 

We start with a pedagogical discussion of X-ray emission from galaxy clusters in order to elucidate important details 
relevant to this work and differences within the literature. We denote the X-ray emissivity [erg s~^ cm~'^ keV~^] at 
the source as 

d^Lx 

'^(-^sourcc; -^sourcc5 '^sourcc)^e; (1) 



dEsQurcc dVs 



where Lx is X-ray luminosity, E is photon energy, V is volume, T = Te is the cluster temperature, Z is metallicity 
relative to solar, Uf. is the electron density, and we use the index "source" for clarity to indicate variables defined 
at the source position (with redshift z). The function A(i?souice; ^source, ^source) [erg s""'^ cm^ keV~^] depends on the 
temperature Tgourcc and the metal abundance ^source of the source, and is given by the theoretical model. 

Consider X-ray photons isotropically emitted from the source of physical size AxAyAi (with i chosen to be along 
the line of sight of the observer at z = 0) with energy range -Esourcc ~ ^'source + Ai?sourcc for a time interval Atsouico- 
Then the total number of the photons is given by 

A/T""' = p ^ ,p "^'^5. X AS,ourceAt,ourccAa;AyAl (2) 

-L^sourcc 0-LJsouice(^ ^^sourcc 

The ideal observer (with perfect efficiency) receives the same number of photons from the solid angle of the source 
Afiobs = AOxAOy = AxAy/d\ for a time interval Afobs with energy range £^obs ^ Eohs + Ai?obs, where those variables 
with index "obs" are defined at the observer's frame, and cJa is the angular diameter distance to the source. If we 
denote the "differential" surface brightness of the source as 

^^theory 

— [erg s"^ cm"^ arcmin"^ keV"^], (3) 

then the number of photons in the observer's frame is given by 

1 dS^^°°^^ 

= 4^dcomi? 1§ AEob^Atob^Aflobs, (4) 
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where c?com is the comoving distance to the source at redshift z. Since the number of photons in both frames is 

y^bs ^ j^soi — — 



invariant, Af^^ = Af]^™'''^, Equations ([2]) and dU are equal and we obtain 



^gthcory ^ ^ ^^^^ d^Lx A^source A^source Aa: Ay 

^-^obs ^^'^com -^source '^-^source'^ Source ^^obs ^-^obs ^^obs 

1 1 d'^Lx 1 



^ (1 + z) dlAe 



47rd2^^ 1+ Z d£'sourccC?V;ourco 1 + Z 

1 fd.V d^L, 



47r(l ^2^) \ ^^com / ^-^sourcc^^ourcc 



1 d'^L, 



47r(l + z)3 dEsomccdVsc 



(5) 



where we use i^souicc = (1 + z)Eohs, Atgourcc = A<obs/(l + z), and d/^ = dcom/(l + z). 

If the profile along the line of sight is properly taken into account, A£ in the above equation will be replaced by the 
integration along the line of sight: 



dEohs 47r(l + z)3 J dEsouiccdVs 
1 



47r(l + z)3 
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^^sourcc^ *^sourcc 

d£ A(-£/sourcc7 -^sourcci -^sourccj^e ' (^) 



Then the surface brightness of the source, 5'^'^°^^ [erg s ^ cm ^ arcmin defined at the observer's frame is given 
by 



theory 



othoory^ f ,p dS^ _ 1 I I a-l.x 

-J dE,^^ ^^^^^ - ^^^^ ^ J di^obs J 



source'^ *^ source 
2 



^ ^ ^-j^ d£ j dEgQ^YCc '^(-^sourcci ^source; ^source)^e 

4 dM*h<=-y(T„)n2, (7) 



47r(l + z) 
1 



47r(l + zY 

where the cooling function A*^'^°''y(Tsourcc) [erg s~^ cm^] is given by 

A '^(^sourcc) — ^-^sourcc ^(-^sourcc ; '^source; '^source) • (8) 

Note that A*^'^°'^y(Tsourcc) indeed depends on ^source as well, but we do not write it explicitly for simplicity of notation. 

2.2. From observed photon counts to surface brightness 

In real observations, we have to take account of the overall response function of the detector and the telescope, 
whose effect is expressed by the energy-dependent effective area A{Eohs) [cm^]. Then the observed number of photons 
per unit time per unit solid angle counted on the detector, d'^M^^ / (dtohsd^obs) [counts s^^ arcmin^^], is modeled as 

^ ■ = / dEobs 1= t| A(Kbs) 



(^^obs^^obs J Eohs dEo\ys 

^ 7 otheory 

— ^^^A(£;obs) 



^source '^^obs 

1 



47r(l + z)3 



yj '^(-^source: ^source? ■^source) . ( -^source \ 
a-Esourco A — . (9) 

-C/source \ + ^ / 

Spectral analysis of the observed X-ray photons yields Tgourcc and ^source that best-fit the data. 
Let us define the effective cooling function as 

V sourccy — 



A{Eohs — Efid) 



/■-E2,ob3(l + z) \{p rp y \ /fP \ 

I 7 7-^ source 7 source^ -^source y . / -'^source \ 
X / dEsource — -A — , (10) 

JEi^absil+z) -C'sourcc \ 1 + Z / 

where i?fid is some fiducial energy, Efi^ = 1 keV for the Chandra analyses presented and compared in this work, and 
i?i_obs and i?2,obs correspond to the observed energy band of the detector (or of the analysis). 
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Before proceeding further, a few comments should be added here. First, this quantity has units of [erg s ^ cm'^ 
keV^^] — [counts s^^ cm^], different from that of A*^''"''^ [erg s~^ cm^]. In addition to covering different energy ranges, 
ROSAT wA Chandra treat exposure maps differently, which leads to different treatments of the "cooling function". 
Therefore one must be careful when comparing emissivity and surface brightness results even in cgs units. Of course, a 
"count" in different observatories corresponds to different energy photons, on average, so "detector" unit results most 
often can not be compared directly either. Choosing Ea^ = 1 keV in BOG is effectively arbitrary but does roughly 
correspond to where Chandra is most sensitive and corresponds to the energy at which the exposure maps used in 
their analysis were computed. We follow suit and adopt i?fid = 1 keV for a more straight-forward comparison. 

Finally the effective cooling function depends on ^source and z in addition to Tgource, but we write it as A°^(Tsource) 
for simplicity. Equations ([9]) and ()10p arc now combined to give 

j2 \f obs 1 1 r 

- / n^d^ X A<=«(T„). (11) 



dtohsdVlohs ^(-Bobs = -Efid) 47r(l + z)'^ 
Comparing this equation with Equation ([7]), we find that if we define the quantity 



dW?-'"" 1 



5*^ = ^ ^ (12) 



dtobsrff^obs ^(-Bobs — -Bfid) 

then we obtain 

= 4.(1 Vz)3 / "'^^ ^ A^^(^"). (13) 

Again it should be noted that S'^^ [erg s""'^ cm~^ arcmin^^ keV^^] {— [counts s~^ cm~^ arcmin"^]) is not directly 

related to the surface brightness 5'^'°°'^, but this notation follows that used in the literature.^ 

In any case, if the effective cooling function is independent of the line-of-sight integral of Equation (fT3|) . such as in 
the isothermal case considered here, we can estimate the density squared integrated along the line of sight as 

1 /■ , S'^'^ [counts s^^ cm^^ arcmin"^! , , 



47r(l + z)4 7 ^ (1 + z)A°ff(T,ourcc) [counts s" 

Even though Sf"" and A°*^(rsourcc) are not the direct observational counterparts of 5^°°''' and A*'^''°''y(rsourco), re- 
spectively, their ratio may be used to estimate the left-hand-side in the above equation. 

2.3. Breakdown of factors important for determining Hq with different calibrations 
As discussed previously, the X-ray emission is 

S'^' ^ j nl K"^ de oc nl A"*^ Ai. (15) 

Since the observed SZE temperature decrement is proportional to the pressure integrated along the line-of-sight, 

ATszE CX /(i/,T,o„co) J "eTsourccrf^ ~ /(!y,T,o„,„,)n-ersourco A^, (16) 

the angular diameter distance, dA, can be estimated by taking advantage of the different dependencies on n^. The 
frequency depen dence of the SZE, f(,^^T^), also depends on temperature when relativistic effects are considered (e.g., 
Iltoh et al.iri998HChallinor fc Lasenbvl 119981 ). Eliminating in favor of Ai yields 

ATszE \'A-ff 



d^cxA^oc — ^. (17) 

V ^source /(i^,T,o„,„) / ^X 

Therefore the estimated Hubble constant should be proportional to 

rp2 cobs j-2 rp2 f2 

TT A/7-1 source "^X ■^(i/^T.oui-oe) ^ source.; (t/,T,„„,,.) 

A<=ff(r_) Aoff(T„)A(£;obs = i?fid)' ^ ^ 

where we assume that the values of ATsze and Af^^ are not affected by the change of the cal ibration. Details on the 
calculation of dis t ances from the an alysis of X-ray and SZE data can be found elsewhere (e.g.. iBirkinshaw et al.llT99TI : 
lBirkinshawlll999HReese et al.ll2002[ ). The dependence of /(i/,Tsourco) is weak and we neglect it in the rest of this 

section, although we do include it in the distance and Hubble constant calculations from the data (Section 

® In R02, some effort was made to distinguish between the observed, denoted as "detector" , and theoretical, denoted as "cgs" , quantities 
and both were presented while the B06 analysis presents only observed (detector) values. 
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Thus the ratio of the estimates of Hq from two different calibrations, 1 and 2, is finally written as 

Hpg^ T| Af (ri)Ai(^fid) 

Ho,i Af{T2)A2{Efid) T2 
Af{Ti)Ai{Efid) 



Tj Af{T2)A2{Efi^)- ^^^^ 
Equation (jl9l) can be rewritten in a form useful in understanding the error budget of the estimate of Hq 

Hog (T2V Kf{T,) AfjT,) A,{En,) 

Ho.i \tJ Af(T2)Af(Ti)^2(i?fid)' ^ ' 

where A2^(Ti) is the effective cooling function using the temperature from observation 1 but the effective area cor- 
responding to observation 2. Because it enters the Hq calculation squared, the effects of the calibration on Te will 
have the greatest impact. For instance, a ^ 10% change in Te would result in a ^ 20% change in the derived Hubble 
parameter. 

The product of the third and fourth factors indicates the effect of the different effective areas evaluated at the same 
temperature (and same abundance). More specifically it is given by 

Af{Ti)Ai{Enci) 
Af{Ti)A2{Efid) 

E2,abs{l+z) A(-Esourcc; Ti,Zi) 

-El,obB(l-|-z) -C/sourcc 



-E2,ob3(l-l-z) A(£'sourcG, T^l, -Z^l) 



7J71 ''\"&UUiHJl^J.5"-L/ A 



Bi.ob=(l + z) 



1 + Z 



(21) 



The above factors may change the inferred Hubble parameter by several percent (see Figure [T] in Section l3.2p . 

3. SPECTROSCOPIC TEMPERATURE OF SZE CLUSTERS WITH Chandra 

We determine the spectroscopic temperature of all 38 clusters used in the BOG analysis with three different versions of 
the calibration database and corresponding ciao versions: calibration versions 3.1.0, 4.1.4, and 4.2.2 with ciao versions 
3.4, 4.1.2, and 4.2. These calibrations will be referred to as 3.1, 4.1, and 4.2 hereafter. 

The most significant change between the 3.1 and 4.1 calibrations is the updated high- resolution mirror assembly 
(HRMA) effective area in version 4.1. A new effective area model was developed after revisiting the ground calibration 
data and ray-trace model, prompted, in large part, by differences in inferred temperatures of massive galaxy clusters 
between Chandra and XMM-Newton. The newer effective area at low energy {E < 5 keV) is < 10% lower than that 
of the older calibration.^ The main change in calibration version 4.2 is the AXAF CCD imaging spectrometer (ACIS) 
contamination model. ^ Based on external calibration source (ECS) measurements, separate ACIS-I and ACIS-S models 
were developed. An offset was added to normal exponential evolution of the contamination model for ACIS-S,"^ while 
there is no such offset applied for ACIS-I. We note that there are a lot of other minor changes between 4.1 and earlier 
versions, and between 4.1 and 4.2. These more minor changes also have the potential to affect changes in the estimated 
cluster temperatures but to a much lesser extent than the changes outlined above. 



3.1. Data Reduction 

All available archival Chandra data for the 38 galaxy clusters in our sample are used in the analysis. This includes 
all data used in BOG as well as observations available at the time but not included and new observations since that 
time that are in the public archive. There are 87 observations among the 38 galaxy clusters, not including two 
observations of A02G7 that have < 1 ks of observatio n time that were not used. The data are summarized in Table [31 
Cluster position, redshift, H I column density (Dic kev fc Lockmanlll990l) . observation identification number, ACIS 
configuration (ACIS-I or ACIS-S), and livetime (effective integration time) are shown. Redshift references may be 
found in BOG. 

The data are reduced with each of the three different versions of the calibration database 3.4, 4.1, and 4.2. The data 
are processed starting with the level 1 events data, removing the cosmic ray afterglow correction, and generating a 
new bad pixel file that accounts for hot pixels and cosmic ray afterglows. Using the newly generated bad pixel file, the 
charge transfer inefficiency correction, time-dependent gain adjustment, and other standard corrections are applied 
to the data. The data are filtered for ASCA grades 0, 2, 3, 4, G and status=0 events and the good time interval 
data provided with the observations are applied. Periods of high background count rate are excised using an iterative 
procedure involving creating light curves in background regions with 259 s bins (following the ACIS "Blank-Sky" 

^ http: / /cxc. harvard.edu/ciao/releasenotes/ciao_4. l_release.html 
^ http: / / cxc.harvard.edu / ciao / releasenotes / ciao_4. 2_release.html 
^ http;/ /cxc. harvard.edu/cal/memos/contam_memo.pdf 
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Background File reduction), and excising time intervals that are in excess of 3 cr (=rms) from the mean background 
count rate. This sigma clipping procedure is iterated until all remaining data lie within 3 a of the mean. The final 
events list is limited to energies 0.7-7.0 keV to exclude the low- and high-energy data that are more strongly affected 
by calibration uncertainties. 

Lightcurve filtering is performed on background regions chosen for each observation. Three circular regions, masking 
out any intervening point sources, are used as the background regions. Investigation of d eep blank-sky exp osures 
shows that front-illuminated (FI) and back-illuminated (BI) chips have different responses (jBonamente et all [2004'). 
In particular, the BI chips are basically constant over the chip but there is a gradient in the FI chips that depends on 
the distance from the readout nodes. For observations taken in the ACIS-S configuration, the cluster observation falls 
on a BI chip. In these cases, background regions are chosen at the periphery of the BI chip, away from the cluster 
emission. For observations taken in the ACIS-I configuration, there are four main FI chips. The cluster falls on one of 
those chips (13). Background regions are chosen from the three remaining chips (10, II, and 12), one region for each 
chip, using locations at the same distance from the readout nodes as the cluster. 

To facilitate point source detection, images and exposure maps are constructed. Images are created by binning the 
data by a factor of 4, resulting in a pixel size of 1".97. Exposure maps are constructed for each observation at an 
energy of 1 keV. For clusters with multiple data sets, the observation with the longest exposure time is used for point 
source detection. A wavelet based source detector is used to find and generate a list of potential point sources. The 
list is examined by eye, removing bogus or suspect detections, and then used as the basis for our point source mask. 

Spectra are extracted from the cluster using a region that encompasses 95% of the galaxy cluster counts, accounting 
for point sources that fall within the region of interest. The same wavelet based source detection algorithm as for 
point source detection is used to find the cluster by having it search for large scales (compared to the point-spread 
function (PSF)) and determine the initial extraction region. The center of the ellipse returned by the source detector 
is adopted as the center of the extraction region and the semimajor and semiminor axes of the ellipse are combined in 
quadrature and its square root used as the radius for the circular extraction region. Background subtracted counts are 
then computed in annuli out to the initial extraction radius and the region encompassing 95% of the cluster counts is 
adopted as the final extraction radius. This provides a formulaic method of constructing spectral extraction regions. 

Calibration version 4.1 is used to construct background and spectral extraction regions that are then used for the 
other two reductions. This provides a uniform reduction and systematic analysis procedure that isolated the effects 
of the calibration. Using these regions, spectra are extracted and responses computed for each observation. Multiple 
observations for a given cluster each have spectra and response files that are then simultaneously fit to a thermal 
spectrum. 



3.2. Spectral Analysis 

Following iBonamente et all ([2006I . [200l . XSPEC (lArnaudI [l99l ID orman fc ArnaudI l2001f ) is used to model the 
intracluster medium (ICM) with a Raymond — Smith spectrum ( Ravmond fc Smith 1977* ) accounting for galactic ex - 
tinction, with solar abundances o f iFeldmanI (.1992) , and cross sections of Il3alucinska-Church fc McCammonI (|1992[ ) 
with an updated He cross section (|Yan et al.lll998D . The analysis uses data in the 0.7-7.0 keV energy range. We use 
the "cstat" statistic in XSPEC, essentially the so-called Cash-statistic ()Cashlll9T9t ). to properly account for low count 
spectral bins. Using the statistic in such cases will result in biased spectral parameters. We perform a Markov 
Chain Monte Carlo (MCMC) analysis, allowing us to compute proper uncertainties on the X-ray cooling function, 
Aeff, which also enters the distance calculation and is derived from the spectral parameters from the MCMC chain. 

A simultaneous fit is performed to all data sets of a given cluster. Responses for each individual cluster observation 
may be quite different even for the same cluster. In particular, the FI and BI ACIS chips show considerably different 
responses. Therefore, for clusters with multiple observations, each individual observation is modeled independently, 
using their individual response and background spectra. Normalizations of the model are allowed to vary for each data 
set while the temperatu re and abundance ar e linke d between the different data sets. Redshift and column density are 
held fixed, adopting the lDickev fc LockmanI (|199O0 Nn values (See Table [3]). 

Markov chains are run for 100,000 iterations. We drop the initial 5000 iterations for the burn-in period but the 
results are insensitive to that choice. The X-ray cooling function, Aoff, is computed at each step in the Markov 
chain, enabling the folding of uncertainties of the spectral parameters to the computed Aeff. Best-fit parameters and 
confidence intervals are computed from the cumulative distribution with 50%, 16%, and 84% probability, corresponding 
to the median and 68% confidence interval. The resultant probability distribution functions for each fit and derived 
parameter of e ach cluster are visually inspected. Convergence and mixing are checked with the Geweke Z-statistic 
(IGewekel[l992h . 

For clusters with multiple observations, a weighted average effective area as a function of energy is used to compute 
the X-ray cooling function. The weighted average effective area, A^^^giE), is the average effective area weighted by the 
livetime (effective integration time) of each observation, Aavg(-E) = ij^^ii-^) * ti)/ where Ai{E) is the effective 
area of observation i, U is the livetime of observation i, and the sum is over the number of observations for that cluster. 
The weighted average effective area is used to compute Aeff at each step in the Markov chain and determine the best-fit 
value and 68% confidence interval. 



3.3. Spectral Results 
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Fig. 1. — Comparison of the effective areas for the 3.1 (red), 4.1 (blue), and 4.2 (black) calibrations for examples of AClS-1 (CL0016+16; 
top) and ACIS-S (MS1054. 5-0321; bottom) observations. The fractional residuals from the 4.2 effective area are shown in the right hand 
panels, AA/A^'^ = (A^ — yl*'^)/yl*'^, where x is either 3.1 (red) or 4.1 (blue). In addition to the changes seen between 1 and 2 keV in 
both ACIS-I and ACIS-S observations, ACIS-S observations also exhibit appreciable differences between the 3.1 and 4.1/4.2 calibrations in 
the 2—5 keV range. The vertical dotted lines mark the 0.7-7.0 keV energy range used in this analysis. 



The newest calibration version, 4.2, will be used as the baseline for comparison with the other two calibrations 
throughout the analysis. We define the mean ratio of parameter P between calibrations as 

/ \ 1 p-^ 

\P^/ = N~i^Pp ^^^^ 

where i denotes the value for each individual cluster, x refers to 3.1 or 4.1, and A'ci is the number of clusters used in 
the calculation. 

Examples of the differences in the effective area for the different calibrations are shown in Figure [T] for both ACIS-I 
(CL0016-1-1609; top) and ACIS-S (MS1054. 5-0321; bottom) observations of galaxy clusters. Shown are the effective 
areas (left) and fractional residuals from the 4.2 calibration (right) for the 3.1 (red), 4.1 (blue), and 4.2 (black) 
calibration versions. The vertical dotted lines denote the 0.7-7.0 keV energy range used in this analysis. Residuals 
are defined as IS.A/A^-'^ = [A^ — A'^-'^)/A'^-'^, where x refers to 3.1 or 4.1. The same definition of residual applies to 
the other cluster properties with A replaced by the parameter of interest. The resulting temperatures follow the same 
trends as the normalization of the effective areas (Tg oc A{E)). 

The results from the Markov chain analysis are summarized in Table |4] for all three calibration versions. The basic 
trend is as expected from the calibration notes and the effective area curves, T^'^ > Tg > T^'^. A2163 appears in the 
table twice. Following the methodology outlined above results in what seems like unrealistically high temperatures, 
16 — 20 keV between the calibrations. In addition, the combination of its formally small uncertainties and it being 
very far from the mean inferred Hubble parameter (due to its high Tg) causes undue infiuence on the resultant Hubble 
parameter (see Section |4]) . This cluster is among a handful of clusters that are known to have H I column densities 
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Fig. 2. — Comparison of different temperature measurements (left) and the corresponding fractional residuals from the 4.2 calibration 
results (right). Temperature determinations using the 3.1 calibration (red) and 4.1 (blue) calibration are compared against temperatures 
using the newest 4.2 version. Also plotted are the ASCA (green) temperatures adopted in R02 for comparison of the 17 clusters that 
overlap with this work. The dotted black line shows the equality relation. Residuals are defined as ATe/T*'^ = (T^ — T^'-^)/T^'^ , where 
X is the 3.1 (red), 4.1 (blue), or ASCA (green) results. Error bars show 68% confidence statistical uncertainties and the A2163 results are 
denoted by stars. Points are slightly offset along the x-axis in both cases in order to make it easier to distinguish between the different 
calibrations. 



that are significantly different from the iDickev fc LockmanI ( 19901 ) values (e.g.. iGovoni et al1l2004[ ). When adopting 
the updated column density of 18.7 x 10^° cm~^ (jGovoni et al.ll2004l) instead of the IDickev k LockmanI ()199Clf ) value 
of 12.1 X 10^" cm~^ the temperatures and other parameters change significantly because the column density strongly 
affects the lower energy part of the spectrum. The updated temperatures are in the range 10 — 14 keV, placing A2163 
more in line with previous results and closer to the mean Hubble constant from the other clusters, lessening the impact 
from this single cli ister. As an aside, a lso note that a newer H I column density study yields an even lower value of 
10.9 X 10^° cm~^ (jKalberla et al.l[2005l ). Choosing an updated H 1 column density destroys our uniform, systematic 
study by treating one clust er special. It is teinpting to t hrow out this cluster to keep the uniformity. We therefore 
present results for both the' Dickev fc LockmanI (|1990( ) and IGovoni et al] (|2004l ) A'h values and also explore the effects 
of removing A2163 from the sample. Initially A2163 has a large impact on the final Hq value. However, once the 
new value of A'h is adopted, A2163 basically has no effect on the results. Table S] shows the results for both column 
density values. The figures show results for only the adopted updated Nn value and a star is used to denote A2163 to 
distinguish it from the other clusters. 

R02 adopted Nn values from spectral fits where available (five clusters), use d the value f rom the Bell Labs H I 
survey (Stark ct al. 1 99^ as adopte d in a detailed analysis for that one cluster (jPonahue et a l. 1999), and adopted 
the DL values (jPickev fc Lockmanlll990t) for the remaining clusters. In particular, the adopted value for A2163 is 
A'h = 16.5 X 10^° cm~^ (jElbaz et al.lll995l ). A quick spectral fit to the Chandra data including A'h as a free parameter 
suggests A'h 16 X 10^° cm~^. Further investigation is beyond the scope of this work. We are concerned with the 
overall differences between calibrations more so than the values of the individual derived quantities for a particular 
cluster. 

Figure [2] shows the comparison of the temperatures from the version 3.1 (red) and 4.1 (blue) calibration versus that 
of the 4.2 calibration (left) along with the fractional residuals compared to the 4.2 calibration results (right). Also 
plotted are the ASCA results (green) for the 17 overlapping clusters used for SZE/X-ray distances using the same SZE 
data as B06 but using ROSAT and ASCA data (R02). The dotted line shows the one-to-one correspondence. There 
is a clear division between the 3.1 and 4.1 results on either side of the equality line, with 3.1 falling above the relation 
and 4.1 falling below, clearly seen in the residuals. A2163 results are shown with stars. 

The mean ratios between temperatures are 



y3.: 



1.06 ±0.05, 



t: 



■4.1 



0.93 ±0.03, 



(23) 



where the above uncertainty is simply the rms in the ratio. The changes between the two older calibrations and the 
newest are roughly the same order 6%) but in different directions. Therefore, the average change in temperatures 
between the 3.1 and 4.1 calibrations is of order ^ 12%, which would produce a ~ 24% difference in the inferred Hubble 
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Fig. 3. — Comparison of different abundance measurements (left) and the corresponding fractional residuals from the 4.2 calibration 
results (right). Abundance determinations using the 3.1 calibration (red) and 4.1 (blue) calibration are compared against metallicities 
using the newest 4.2 version. Also plotted are the ASCA (green) abundances adopted in R02 for comparison of the 17 clusters that overlap 
with this work. The dotted black line shows the equality relation. Error bars show 68% confidence statistical uncertainties and the A2163 
results are denoted by stars. Points are slightly offset along the x-axis in both cases in order to make it easier to distinguish between the 
different calibrations. 
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Fig. 4. — Comparison of different effective X-ray cooling functions (left) and the corresponding fractional residuals (right). Effective 
cooling function determinations using calibration versions 3.1 (red) and 4.1 (blue) are compared against emissivities using the new 4.2 
calibration results. The dotted black line shows the equality relation. Error bars show 68% confidence statistical uncertainties and the 
A2163 results are denoted by stars. Points are slightly offset along the as-axis in both cases in order to make it easier to distinguish between 
the different calibrations. 



The mean ratios comparing against the BOG results are 



/in. 



= 1.05 ±0.11, 



\TB06 



= 0.92 ±0.10, 



rp4 



= 0.99 ±0.11. 



(24) 



Therefore it is likely that there is a small (~ 5%) overall systematic between the spectral analysis of B06 and this 
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TABLE 1 

Compilation of Mean Ratios: Updated A2163 A^h 



Parameter 3.1/4.2 4.1/4.2 3.1/B06 4.1/B06 4.2/B06 A5CA/4.2 

Te 1.06 ±0.05 0.93 ±0.03 1.05 ±0.11 0.92 ±0.10 0.99 ±0.11 0.98 ±0.12 

Z 1.08 ±0.21 0.96 ±0.04 1.16 ±0.43 1.03 ± 0.31 1.08 ± 0.34 0.66 ± 0.28 

Aoff 1.01 ±0.01 1.01 ±0.01 1.03 ±0.07 1.03 ± 0.07 1.02 ± 0.07 

/(i.,Te) 0.998 ±0.002 1.002 ±0.001 0.999 ± 0.003 1.003 ± 0.004 1.001 ± 0.003 

A(ikev)=' 1.01 ±0.02 0.95 ±0.01 0.96 ±0.10 0.91 ±0.10 0.95 ±0.10 

dj^ 0.93 ±0.08 1.13 ±0.06 1.06 ± 0.24 1.29 ± 0.29 1.15 ±0.25 1.07 ±0.37 

B06 are the effecti ve areas from the 3.1 calibration using only those data sets that appear in 

IBonamente et al.|[2006l . 

B06 are the published distances from IBonamente et al.ir2006l . 



work. Given the various decisions on background regions, spectral extraction regions, and so forth, this is reasonable 
agreement. This analysis removes those uncertainties between analyses by performing the same systematic analysis on 
the same data for three different Chandra calibration versions. We focus on the differences between the calibrations 
rather than the values themselves. The mean ratios of the temperatures as well as other parameters considered here 
are summarized in Table [TJ 

The abundances from the spectral fits (left) and the fractional residuals with the 4.2 results (right) are plotted in 
Figure|3l Again calibration versions 3.1 (red) and 4.1 (blue) are plotted against the 4.2 calibration and the dotted line 
shows equality. Also shown are the ASCA abundances adopted in R02 . There is no clear trend in the Chandra results 
as is the case for the temperatures. There is a mild offset between the Chandra and ASCA metallicities. The ASCA 
results are compiled from the literature and the data are not uniformly analyzed, the abundances used for the analyses 
often differing. The abundances do change the effective cooling functions, Aoff. However, the effects of metallicity on 
Aoff are small, typically on the order of ^ 1% and < 5% even for changing the abundance by factors of 2 or 3. 

Derived effective cooling functions from the spectral analysis (left) and corresponding fractional residuals from the 
4.2 calibration (right) are illustrated in Figure ID Results for calibration versions 3.1 (red) and 4.1 (blue) are shown 
against that of version 4.2. The equality relation is also shown. There is a clear trend for Aes to be greater in 3.1 and 
4.1 compared to the 4.2 results but as shown by the residuals, it is a small effect, < 4% for 3.1 and < 2% for the 4.1 
results. 

The abundance of metals does not directly enter the distance calculation. However, it indirectly enters through the X- 
ray cooling function, Agg. The X-ray temperature also indirectly enters the calculation through relativistic corrections 
to the frequency dependence of the SZE, /(j^.tj, where v is the frequency of the observations (see Section IJ]). Mean 
ratios and r ms's for the temperature, abundance. X-ray cooling function, SZE frequency function including relativistic 
corrections (|Itoh et al...l998.) . the effective area at 1 keV, and angular diameter distance (see Section|4]) are summarized 
in Table m 



4. IMPLICATIONS TO THE ESTIMATE OF Hq 

Direct angular diamet er distances, dA, from a combined SZE and X-ray analysis are stra ight forward to compute, in 
theory, (see, for example iBirkinshaw et al.lll991l : iReese et al.ll2002l : IBonamente et al.l[2006[ ). particularly for the simple 
isothermal /3-model. The actual observations and parameter extraction from the data are not as straight forward. We 
concern ourselves here with only the components of the distance calculation that involve X-ray spectral properties or 
are derived from them. SZE/X-ray derived distances have the following dependence on X-ray spectral properties 

, Aoff ActfA{Efid) 

dA oc — oc -—^ , (25) 

where Tg is the electron temperature, Acg is the effective X-ray cooling function (see Section [2^ . and f( u,T^) is the 
spectral depen dence of the SZE at frequency, including relativistic corrections, which depend on Tg (e.g.. lltoh et al.l 
[1998: Challinor fc Lasenbv l[l998.). Sx is the X-ray surface brightness that depends on the effective area that changes 
between the calibrations, A{E) is the effective area of the observatory, and we have used the fact that Sx oc TVx /A{Efid), 
where J\fx is the number of observed counts and is constant. Since the exposure maps are all computed at 1 keV, Sx 
simply changes by the ratio of effective areas at E^d = 1 keV. It is not immediately obvious how the spectral results 
will affect the final Hubble constant due to the complexity of the dependencies on spectral parameters. However, it is 
clear that a ~ 10% change in will have an appreciable effect on the distances and therefore on the inferred Hubble 
constant. 

Cluster spatial properties from /3-model fits (B06) are adopted for the distance calculation. X-ray spectral properties 
from this analysis for the three calibrations are combined with the adopted cluster spatial properties in order to compute 
distances to each galaxy cluster. We estimate the uncertainty on dA by backing out the X-ray spectral variable 
uncertainty from the published uncertainty and including the new Chandra calibration spectral result uncertainties 
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assuming everything adds in quadrature. Namely we compute the uncertainty in cIa by 

where x refers to one of the three calibrations of this work and 506 refers to the published values (B06). The average 
of the positive and negative uncertainties is used for the uncertainty. This method preserves, as best we can, the 
correlations amon g the parameters. In particular, rc and /3 from the /?-model are strongly correlated (for an example 
in this context see lReese et al.l[200Cl[ ). However, the spectral results are independent of the /3-model. Therefore we can 
preserve this correlation through uncertainty propagation with this method. The temperature uncertainty comes with 
a factor of 2 because the angular diameter distance is inversely proportional to temperature squared, dA oc T~^. We 
include the additional sources of statistical uncertainty, ~ 19%, from Table 3 of B06 by adding in quadrature to the 
uncertainties computed from Equation (j26p . These total statistical uncertainties are then used for the Hq calculation 
and are summarized in Table|4l which shows the derived angular diameter distances with 68% statistical uncertainties. 

The Hubble constant is computed by performing a fit to the cluster distances using the theoretical angular 
diamete r distance relation fo r a flat , A-dominated universe with i^A = 0.73 and = 0.27 consistent with the WMAP 
results (jKomatsu et al.ir2010L 120091 ). The resulting angular diameter distances for the full sample and updated iVn for 
A2163 yield 

i?^-i = 70.0±3.7 km s-i Mpc-i, (27) 
i?(^-^ = 55.4±2.9 km s-i Mpc-i, (28) 
77(^-2 = 63.7 ±3.3 km s-i Mpc-i, (29) 

where the uncertainties are statistical only at 68% confidence with = 40.6, 34.4, 38.8 for 37 degrees of freedom for 
the 3.1, 4.1, and 4.2 calibrations respectively. Solely from the changing calibration, taking the 4.2 calibration as the 
baseline, H^-^/H^-^ = 1.10 and H^-^/H^-^ 0.87, showing changes of roughly 10% and 13% in the determination of 
the Hubble constant, respectively. There is a 23% change in Hq between the 3.1 and 4.1 calibrations, consistent 
with the expectation from the mean of the ratio of temperatures (see Section [331 and Table [T]). 

In Figure [SJ we plot the angular diameter distances (top) and corresponding Hubble constants (bottom) determined 
using the results from the 3.1 (red), 4.1 (blue), and 4.2 (black) calibrations for each cluster. In both cases, the right 
panels show the fractional residuals from the 4.2 calibration results. Results from R02 are also plotted for comparison 
and denoted as ASCA (green). The best- fit theoretical angular diameter distance relations (top) and corresponding 
Hubble constants (bottom) are also shown. Angular diameter distances are also summarized in Table |4] for each of 
the calibrations. Uncertainties are 68% confidence and include both the uncertainty from the fit and the additional 
sources of statistical uncertainty from Table 3 of B06 . In general, the distances follow the trends expected from the 
temperatures, namely, < d\'^ < d\^, indicating that the temperature changes dominate the changes in distances. 

Hubble parameters are computed from the determined angular diameter distances using for the original Nn = 
12.1 X lO^'' cm^^ for A2163, the updated A2163 iVn, and excluding A2163 from the fit for each of the three calibrations. 
The results are compiled in Table[2l where the uncertainties are 68% confidence statistical uncertainties and the number 
in parentheses is the the best-fit value. Hubble parameters are computed from fits for the 17 clusters that 

overlap in the R02 sample, becoming 16 clusters when A2163 is excluded. A simple average is also computed for each 
case by first converting the c?a for each cluster into a Hubble parameter and then averaging. Finally, a fit is done 
for the published values of B06 and R02 both including and excluding A2163. 

Table [2] summarizes the full set of Hubble parameter estimates. The beginning of each section marks the Hq 
estimation method (x^ or average) followed by the number of clusters used in each calculation. The last two sections 
are the results of performing x^ fits both including and not including A2163 using the published distances from B06 
and R02. In the case of R02, the first entry includes the entire sample of 18 clusters although only 17 of which overlap 
with the B06 sample. A fit of just those 17 appears in the second column and then the results without A2163 in the 
third column. 

Using the means of the ratios of the X-ray spectral parameters with respect to the calibration version 4.2 results 
(Table [J) and scaling H^ '^ implies H^-^ = 69.2 km s"^ Mpc"^ and H§-^ = 56.5 km s"^ Mpc~\ very close to the 
Hubble p arameters from the ana lysis. This is also true when dropping A2163 from the analysis but not true when 
using the lDickev fc LockmanI (|1990r ) A^h value for A2163. In that case, using the ratios to scale the 4.2 results predicts 
ijs i = 75.5 km s"^ Mpc"i and i?^-^ = 62.3 km s^^ Mpc^^ compared to the results 82.8 and 58.4 km s"^ Mpc~\ 
respectively. 

5. DISCUSSION 

We perform a uniform, systematic spectral analysis of 38 galaxy clusters using three different Chandra calibrations 
and find significant differences in the inferred spectral properties of galaxy clusters between the calibrations. Using 
the newest calibration, 4.2, as the baseline to which to compare, the temperatures change ~ 6% on average in the 
3.1 and 4.1 calibrations, but in opposite directions. In particular, the temperature changes between the extreme cases 
(3.1 and 4.1) show a ~ 13% difference in Tg, on average. These results are consistent with an analysis of 10 galaxy 
clusters that found the 4.1 calibration results yield ~ 10% lower temperatures than the 3.1 calibration (|Gaeta,2009.) . 
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Fig. 5. — Angular diameter distances (top) and corresponding Hubble parameters (bottom) determined using the results from the 3.1 
(red), 4.1 (blue), and 4.2 (black) calibrations. In both cases, the right panels show the fractional residuals from the 4.2 calibration results. 
Results from R02 are also plotted for comparison and denoted as ASCA (green). The bestOfit theoretical angular diameter distance relations 
are also shown for each calibration (top, left) with corresponding Hubble parameters in the legend. The best-fit Hubble parameters are 
shown in the Hg figure (bottom, left) and corresponding residuals (bottom, right). Error bars show 68% confidence statistical uncertainties 
and the A2163 results are denoted by stars. Points are offset slightly along the x-axis in order to make it easier to distinguish between the 
different calibrations. 



The differences in spectral properties of galaxy clusters between the calibrations may have profound implications on 
inferred cosmological parameters from galaxy cluster studies. 

Using the simple isothermal /3- model as a vehicle we explore the ramifications of the effects of the Chandra calibration 
on a particular cosmological application, distances to galaxy clusters from a combined analysis of SZE and X-ray data. 
Because Hq oc T^, changes in cluster temperatures have a potentially large impact on the inferred Hubble parameter. 
The 3.1, 4.1, and 4.2 calibrations imply Hubble constants of 70.0, 55.4, and 63.7 km s~^ Mpc~^ using the updated 
column density for A2163 results. The results remain essentially unchanged when removing A2163 from the sample to 
preserve the uniformity of the analysis to the last detail. Simply from the effects of the Chandra calibration, the 3.1 
results yield a 10% higher Hq than the most recent 4.2 calibration and the 4.1 results yield a 13% lower Hq than the 
4.2 calibration results. This is in rough agreement with what the mean ratios of temperatures would predict, 13% and 
14% effects for the 3.1 and 4.1 calibrations, respectively. This strongly suggests that, although there are a number 
of quantities that change with the spectral results, the change in Te is, by far, the most important when considering 
SZE/X-ray derived distance based estimates of the Hubble parameter. In the most extreme case, comparing the 3.1 
to the 4.1 results, there is a ~ 24% change in the Hubble parameter due to the change in Chandra calibration. 

Although the isothermal assumption is over-simplistic, it is sufficient to study the effects of the new calibration on 
SZE/X-ray derived distances. B06 showed that the Hubble parameter estimates from SZE/X-ray distances do not 
depend as strongly on the model as one might naively believe. In particular, the isothermal /3-model yields results 
consistent with the more sophisticated hydrostatic equilibrium model and an isothermal /3-model with the central 
regions of the X-ray data removed that they considered. This simple isothermal /3-model also facilitates isolating the 
effects that the different spectral parameters have on the distances and therefore on the Hubble constant. 

SZE/X-ray distances tend to favor a Hubble constant of order 60 km s~^ Mpc~^ (e.g., .Carlstrom et al.l [20021 : 
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TABLE 2 



Compilation of Hq Results 




DL A2163 A^H 


Updated A2163 A^h 


No A2163 


full sample 

rr4.2 

Jh 


38 

82.8 ± 4.6(75.9) 
58.4 ±3.1(42.4) 
68.8 ±3.7(52.2) 


38 

70.0 ± 3.7(40.6) 
55.4 ±2.9(34.3) 
63.7 ±3.3(38.8) 


37 

69.7 ± 3.7(40.5) 
55.5 ± 2.9(34.3) 
63.7 ±3.4(38.8) 


X R02 overlap 


17 

90.1 ± 7.0(48.5) 

58.2 ±4.1(21.2) 

'zn 1 -1- c; 1 /OQ '7^ 
/U.i it 0.1(ZO. ( ) 


17 

66.9 ± 4.7(15.8) 
52.3 ±3.6(11.8) 

tin ^ J- A Q/1/1 /l^ 
DU.O ± 4.o(14.4J 


16 

66.1 ± 4.8(15.5) 

52.2 ± 3.8(11.8) 

oU.z iL 4. 4(^14. 


Avg lull sample 

rr3.1 


38 

66.1 ±30.3 
52.5 ± 17.2 
59.7 ±22.1 


38 

62.9 ±21.7 
51.3 ± 15.4 
58.0 ± 18.9 


37 

62.6 ±21.9 
51.2 ± 15.6 
57.8 ± 19.1 


Avg K,02 overlap 

rr4.2 


17 

68.3 ±36.1 
51.8 ± 16.9 
60.0 ±23.3 


17 

61.2 ± 17.6 
49.2 ± 12.0 
56.1 ± 15.5 


16 

60.4 ± 17.7 
48.9 ± 12.3 
55.6 ± 15.9 


B06 refit 
-Jlil 


38 

76.2 ±4.1(55.9) 




37 

73.5 ±4.1(51.7) 


X^ R02 refit 


18 

60.8 ±4.0(16.5) 


17 

60.5 ±4.1(16.4) 


16 

60.7 ±4.3(16.4) 



Note. — DL is the lDickev fc Lockmaiil | |T99CD H I survey. The first row of 
each section describes the type of Hq calculation (x'^ or average) followed by 
the number of clusters used in the fit for each of the three cases considered. 
The numbers in parentheses are the x^ value at the best fit. The last two 
sections are fits to the published B06 and R02 distances. For R02, the first 
column is the full sample of R02, the second column include the 17 overlapping 
clusters with this work, and the third column then also excludes A2163 from 
the x^ fit- 



iReese et al.ll2002l ) with a few exceptions (e.g.. iBonamente et al.|[2006l : iMason et al.ll2001| ). This is in contrast to other 
probes that f avor Hp ~ 70 km s ~^ Mpc~^, such as the HST Hp key project ( Frccdman ct al. 2001), recent supernova 
results (e.g.. iHicken et al1l2009[ ). and recent WM^AP results (jKomatsu et al. 20 0ft. i201Qi). Potential systematics in 
SZE/X-ray derived estimates of Hp are still for midable (see for example, Birkinsh aw et al.lll991l : lHughes fc Birkinshawl 
ll998tlReese et all2002l : IBonamente et al.ll2006[ ). making the results consistent within the uncertainties. It is still curious 
that SZE/X-ray distances tend to favor a lower Hubble constant. 

One possible explanation of this bias in Hp from SZE/X-ray determined distances to g alaxy clusters is the corn bined 
effects from the degree of inhomogeneity and the multi-temperature structure of the ICM (jKawahara et al.ll2008al ). The 
bias vanishes in the limit of an isothermal and homogeneous ICM, an idealized and unrealistic limit. The presence of 
inhomogeneity in the cluster gas is suggested by a detailed study of the nearby galaxy cluster A3667 (jKawahara et al.l 
I20081i ). which finds 30%-40% rms density fluctuatio ns in that cluster. In addition, 20% rms temperature fluctuations 
have been seen in Hydra A ('Si mionescu et alT 2009V The theoretical underpinning of the inhomogeneity model was 
boosted because the nearby galaxy cluster A3667 exhibits the expected lognormal signature (Kawahara et al.i, 2008b) . 

This model implies that the bias in the Hubble constant can be decomposed into three factors 

r _ -ffo,cst /rp \XT{Tspccj /on\ 

where Xa: Xt(7cw), and Xt(7spcc)/xt(7cw) represent the systematic errors due to the presence of ICM gas inho- 
mogeneities, non-isothermality, and the difference between the spectroscopic (Tgpcc), and emission-weighted (Tew) 
temperatures. Numerical values for each of the above bias factors, of course, depend crucially on the degree of inhomo- 
geneities and the temperature structure of the ICM. However, simulated clusters suggest that these biases are roughly 
(10 — 30)% overestimate, (0 — 20)% underestimate, and (10 — 20)% underestimate, respectively, resul ting in an overall 
(10 — 2 0)% underestimate in the Hubble parameter. For analytic expressions and further details, see IKawahara et "all 
(l2008aD . 

We note that the above result is consistent with the many studies of the statistical and systematic uncertain- 
ties in Hubble constant determinations from SZE/X-ray dis tance s (e.g.. | Inagaki et al. 1995i: IKobavashi et al.lll996l : 
lYoshikawa et al.lll998t IHughes fc Birkinsh aw 1998; Sulkaneiil Il99l IWang fc Fanll2006D that find only smah biases if 
any at all. All of these studies use an emission- weighted temperature for the simulated clusters. This corresponds to 
neglecting the third factor in Equation (1501). and thus leading to no substantial bias because the first and second factors 
coincidentally compensate each other. iMazzotta et al.l (|2004[ ) are the first to point out clearly that the spectroscopic 
temperature, Tgpec, is systematically lower than the emission- weighted temperature. Tew (see also [Mathiesen fc EvrardI 
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l200lURasia et al.ll2005l) . so that the third factor in Equation ^ is essential (|Kawahara et aLll2007l . l2008aD . 

This bias can theoretically be accounted for with more realistic models of the ICM. For example, once the variance 
of the ICM inhomogeneity is known, a fit to even a simplist ic model for the temperature profile largely avoids the 
bias in Hq from SZE/X-ray distances (jKawahara et alJl2008al ). Progress on this front will require a multi- wavelength 
approach combining recent SZE experiments, deep X-ray observations, and weak lensing measurements. 

The large potential systematics in SZE/X-ray derived distances and inferred Hubble constant meai is that the various 
results are consistent w ith each other and other probes of Hq (for detail of the systematics, see, e.g.. lReese et al.ll2Q02l : 
iBonamente et alll2006[ ). Here, we took a detailed look at one of those systematics, the effect of the Chandra calibration 
on cluster temperatures. There have been very little work on including calibration uncertainties in X-ray analysis. 
However, a Monte Carlo approach to incorporate calibration uncertainty for parameter estimation from Chandra 
ACIS-S observations has been developed (Drake ct al. 2006). Better methods of incorporating uncertainties, especially 
systematics such as instrumental calibration, will need to be developed. 

There are alternative methods to X-ray spectroscopy for determining galaxy cluster temperatures. Theoretical studies 
of non-parametric deprojection methods of SZE and X-ray imaging data on both idealized and sim ulated clusters sug- 
gest that cluster tem perature profiles may be accurately re constructed without X-ray spectroscopy ([Yoshikawa &: Sutol 
119991 iPuchwein fc B artclmann 2006; lAmeglio et al.l[2"007| ). Both parametric and non-parametric methods applied to 
actual SZ E and X-ray data also show b road agreem ent with tem perature profiles derived from X-ray spectroscopy 
(iKita vama et all 120041 : iNord et al.l 120091 : [Mroczkowski et all l2009t ) . Cluster tei nperatures may also be inferred from 
SZE data only if one assumes the value for the gas mass fraction of the cluster ([LaRoque et al.ll2Q06l : iJov et al.ll2001|) . 
Because these methods use only X-ray imaging data, they do not require the longer exposure times necessary for spec- 
troscopic Te measurements. In addition, these alternative temperature measurements may, in theory, alleviate some 
of the dependence of derived cosmological parameters on the Chandra calibration. However, most methods currently 
have large uncertainties on the cluster temperature and are complicated by the fact that they still depend on the 
Chandra (or other X-ray observatory) calibration through the effective area for both exposure maps and the cooling 
function calculation. 

The determination of galaxy cluster temperatures is particularly important because it is widely used to infer the 
gravitational mass of clusters. Potential systematics on cluster temperatures will have a strong impact on cosmological 
parameters when using clusters as probes of cosmology. In particular, temperature systematics will strongly affect the 
normalization of the matter p ower spectrum, as, through the cluste r temperature— halo mass relation in conjunction 
with cluster abundances (e.g.. lRasia et al.l[2005HShimizu et al.ll2006( ). 

Analysis of observations of galaxy clusters does have the potential to provide insight into cosm ology and has bee n 
successful in the past, favoring low flM long before hints of a cosmological constant appeared fe.g.. lWhite et al.lll993D . 
Surveys of galaxy clusters have the tantalizing appeal that they probe the growth of structure, one o f the few probes 
to do so, and have the potential to constrain tightly the equation of state of the dark energy fe.g.. iBartlett fc Silkl 
[1994; Holder ct al. 2000; Haiman et al. 2001; Maiumdar & Mohr 2004). However, in order to realize that potential, 
both galaxy clusters themselves and the details of the instruments must be understood precisely which will be feasible 
with combined efforts from ongoing observations and planned missions and observatories covering a wide range of 
wavelengths. 
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37 


47.8 


0.583 


4.96 


551 I 


44 




















1667 I 


44 


MACSJ2129.4-0741 


21 


29 


26.0 


-07 


41 


28.7 


0.570 


4.84 


3199 I 


18 




















3595 I 


19 


RXJ2129.7+0005 


21 


29 


39.9 


+00 


05 


19.8 


0.235 


4.28 


552 I 


10 


MACSJ2214.9-1359 


22 


14 


57.3 


-14 


00 


12.3 


0.483 


3.28 


3259 I 


19 




















5011 I 


18 


MACSJ2228. 5+2036 


22 


28 


33.0 


+20 


37 


14.4 


0.412 


4.58 


3285 I 


20 



^ IDickey h LockmaTI U55^) values. 

^ Both values are considered because A2163 is known to have an iVfj significantly different from rTTT^ey & l.ockman| j 15501 '). 



TABLE 4 
Spectroscopic Results 



00 



(keV) 



2 



3^ 

(Gpc) 



(keV) 



2 



(Gpc) 



(keV) 



2 

(^0) 



(Gpc) 



CL0016 + 1609 


9 


-0 


51 
.47 





.,37ll!! 


.081 
.079 


2. 


.502 + ;!' 


.030 
.028 


1. 


.632 + []' 


.400 
.400 


8. 


,9912 


39 
.36 





34512; 


A0068 


11 


.28i; 


64 
.19 





.4751!! 


.239 
.218 


3. 


.018 + ^' 


.099 
.088 


0. 


.509 + 2 


.203 
.203 


8. 


,5712 


82 
.74 





,508 + 2' 


A0267 


6 


■■"-0 


32 
.30 





.530 + !! 


.105 
.104 


2. 


.976 + ^ 


.050 
.048 


0. 


.978 + 2' 


.303 
.303 


6. 


,2412 


29 
.27 





,482 + 2' 


A0370 


8 


.rs+_» 


42 
.40 





.3751- 


.079 
.075 


2, 


.263"^Q 


.028 
.026 


1, 


.707+2' 


.504 
.504 


7. 


90+0 
'" -0 


37 
.34 





40312' 


MS0451. 6-0305 


10 


.^«l°o 


55 
.53 





.504l0„ 


.094 
.091 


2, 


179+2' 


029 
.028 


1, 


.e4i+g' 


.401 
.401 


9. 


,5912 


50 
.43 





386 + 2' 


MACS,J0647.7+7015 


12 


.07+1- 


22 
.00 





.367+!! 


.144 
.138 


2. 


633+2' 


.053 
.049 


0, 


.858+0. 


.271 
.271 


10. 


,0612 


SO 
.69 





,25912: 


A0586 


7 


16 + 1 


43 
.39 





.659 + !! 


.141 
.127 


2. 


,956 + 2 


.074 
.060 


n. 


.598 + 2' 


.197 
.197 


6. 


.,, + 
■'' -0 


34 
.32 





,58212; 


M ACS J0744. 8+3927 


8 


99 + 


44 
.39 





.4361!! 


.086 
.080 


2, 


,462 + 2 


.035 
.033 


1. 


.667 + 2' 


.485 
.485 


8. 


,0312 


36 
.34 





,379 + 2' 


A0611 


6 


89+"' 


27 
.27 





.393+!!' 


.076 
.073 


2, 


,47812 


.031 
.029 


0. 


.852 + 2' 


.254 
.254 


6. 


,2812 


25 
.24 





'i54+0. 


A0665 


8 


.47t»„ 


27 
.28 





.359 + !!' 


.061 
.061 


2, 


983 + 2 


.027 
.026 


0. 


.845t2' 


.233 
.233 


7. 


30+0 


19 
.18 





,319 + 2' 


A0697 


10 




59 
.57 





.518l» 


.206 
.101 


3. 


,013 + 2 


.044 
.042 


0. 


.872 + 2' 


.266 
.266 


9. 


26 + 


50 
.41 





,44912: 


A0773 


8 


+ 0, 


32 
.32 





.509 + !!' 


.079 
.076 


3. 


14512 


.037 
.036 


1. 


.509 + 2 


.438 
.438 


7. 


,2512 


24 
.21 





,43512; 


ZW3146 


6 




10 
.10 





.477l»„ 


.030 
.029 


2, 


81512 


.014 
.014 


1, 


.048 + 2 


.329 
.329 


6. 


,ool2 


10 
.09 





,44712; 


MS1054. 5-0321 


11 


.301- 


30 
.99 





.134+0. 


.132 
.124 


1, 


862 + 2 


.034 
,027 


1. 


.378 + 2' 


.423 
.423 


9. 


79+0 


84 
.75 





,14212' 


MS1137. 5+6625 


6 


■^*t"o 


64 
.55 





.376+!!' 


.184 
.167 


2, 


,052 + 2 


.062 
.057 


2. 


.586 + 2' 


.919 
.919 


6. 


05 + 


.53 
.45 





'i59 + 0. 

.'■'' -0 


M ACS Jl 149. 5+2223 


9 


.80t!J 


77 

.68 





.244 + !!' 


.117 
.115 


2, 


648 + 2 


.048 
.045 


1, 


.410+2' 


.423 
.423 


8. 


,9212 


63 
.55 





232 + 2' 


A1413 


7 


.65tS 


12 
.11 





.4851S 


.028 
.028 


3. 


295+2' 


.014 
.014 


0. 


.588+2' 


.203 
.203 


6. 


,8612 


08 
.08 





.450+2: 


CLJ1226. 9+3332 


14 


.18+1 


33 
.26 





.11312: 


,139 
.110 


2. 


,33712; 


,046 
.033 


1. 


,007+g: 


,338 
.338 


11, 


.951 + 


.14 
.91 





.169+2' 


M ACS J131 1.0- 0310 


6 


.5o+o: 


26 
.24 





.431+0; 


066 
.065 


2. 


,661+2 


,032 
.032 


2. 


,115+0; 


,792 
.792 


5, 


5512 


.17 
.16 





.45612: 


A1689 


11 




12 
.12 





.324+2; 


,022 
.021 


3. 


,18912: 


,009 
.009 


0. 


.559+2; 


,141 
.141 


9, 


'7612' 


.12 
.11 





.268+2' 


RXJ1347. 5-1145 


14 


.35±g: 


39 
.39 





.553+0 


,049 
.048 


2. 


,598+2 


,017 
.016 


0. 


.870+0. 


,267 
.267 


11. 


.81I2 


.29 
.23 





.457I2: 


MS1358. 4+6245 


7 


■«i°o: 


36 
.29 





.52510. 


.085 
.082 


2, 


433 + 2 


.033 
.032 


1. 


.324 + 2 


.464 
.464 


6. 


57+0 


.25 
.24 





.50112; 


A1S35 


8 




08 
.08 





.455l0„ 


.017 
.016 


2, 


940l2 


.007 
.007 


0. 


.695 + 2 


.168 
.168 


7. 


14+0 


05 
.05 





.44412; 


MACS J1423. 8+2404 


6 


.25l0„ 


15 
.15 





.57810. 


.048 
.045 


2, 


269 + 2 


.022 
,019 


2. 


.243 + 2' 


.841 
.841 


5. 


77+0 


.12 
.11 





.623+2' 


A1914 


11 




36 
.36 





.364+!;' 


.065 
.063 


3. 


,11612 


.027 
.026 


0. 


.596+2' 


.158 
.158 


9 


'^^tS 


.29 
.24 





.372 + 2' 


A1995 


9 


.ooiS 


30 
.30 





.47610. 


.066 
.064 


2, 


606+2' 


.025 
.024 


1, 


034^0. 


.246 
.246 


7. 


,7812 


.25 
.23 





.422 + 2' 


A2111 


7 


.69i!! 


94 
.78 





..364 + !!. 


.219 
.194 


2, 


813+2' 


.097 
.085 


0, 


.903+2' 


.440 
.440 


7. 


31 + 

■"--0 


91 
.73 





.169+2' 


A2163'' 


20 


.98i° 


70 
.69 





.27710 


.065 
.065 


2. 


,64412 


.022 
.023 


0. 


,247 + 2 


.065 
.065 


14. 


•"'-0 


34 
.34 





,34412: 


A2163*^ 


12 


.381!; 


33 
.31 





.319 + !! 


.037 
.037 


2. 


,54512 


.013 
.012 


0. 


,64112 


.168 
.168 


10. 


,2412 


19 
.18 





,318 + 2' 


A2204 


8 


.061!! 


.09 





.57812 


.021 


3. 


,034 + 2 


.010 


1. 


,224 + 2 


.333 


6. 


,4912 


.06 
.05 





,55612: 


A2218 


7 


.23lS 


13 
.18 





.269 + 2 


.056 
.054 


2. 


,87512 


.026 
.024 


1. 


,52412 


.405 
.405 


6. 


'= = tS 


.17 
.17 





263 + 2' 


RX.I1716. 4+6708 


6 


.29lg 


78 
.66 





.867+2' 


.319 
.263 


2, 


21512' 


.127 
.108 


2, 


222II 


.134 
.134 


5. 


9612 


.76 
.62 





,748+2' 


A2259 


5 


.35l0„ 


28 
.27 





.400+2' 


.129 
.116 


2. 


.992+^' 


.084 
.069 


0. 


629 + 2 


.457 
.457 


4. 


.8812 


.26 
.24 





,325 + 2' 


A2261 


8 


.361!! 


28 
.28 





.51412 


.065 
.063 


3. 


.022 + ^- 


.029 
.028 


0. 


765 + 2' 


.227 
.227 


7. 


.3312 


.21 
.19 





.519I2: 


MS2n53. 7-0449 


5 


■ - 


95 
.76 





.2.34 + 2 


.206 
.178 


2. 


.221 + 2 


.106 
.097 


2. 


343+!;' 


.025 
.025 


5. 


,50l2 


.91 
.68 





,180+2' 


MACS,J2129, 4-0741 


9 




79 
.67 





.646 + 2 


.145 
.138 


2. 


714 + "' 


.061 
.056 


1. 


,10512 


.356 
.356 


8. 


.6012 


.62 
.59 





,48712: 


R.X ,72129. 7 + 0005 




.561!! 


23 
.22 





.722 + 2 


.129 
.118 


2. 


.910 + S 


.071 
.067 


1. 


02512' 


.376 
.376 


5 


.3612 


.22 
.22 





,62712; 


MACS,I2214.9-1359 


10 


.181!! 


81 
.73 





.422+2' 


.122 
.123 


2, 


.794+0' 


.050 
.048 


1, 


933+2' 


.541 
.541 


8, 


.1612 


.53 
.49 





.388+0' 


M ACS J2228. 5 + 2036 


8 


■«t°o: 


69 
.63 





.552+2; 


.145 
.134 


2. 


,8i6lg; 


.065 
.060 


2. 


,552+2; 


,762 
.762 


7, 


«t°o 


.60 
.50 





.449+2: 



.067 
.065 
.176 
.160 
.090 
.086 
.080 
.079 
.076 
.074 
.107 
.105 
.120 
.110 
.072 
.068 
.069 
.063 
.051 
.050 
.088 
.083 
.063 
.060 
.027 
.026 
.109 
,105 
.168 
,153 
,094 
,094 
,024 
.023 
.125 
.123 
.064 
.059 
.017 
.016 
.038 
.037 
.072 
.070 
.015 
.014 
.044 
.042 
.052 
.051 
.051 
.049 
.177 
.155 
,040 
.041 
,028 
,027 
,019 
,01s 
,046 
.046 
.278 
.229 
.113 
.104 
.058 
.057 
.210 
.163 
.123 
.114 
.114 
.105 
.092 
.090 
.124 
.115 



2.518 



+ 0. 
-0 



3,066 
2,976 
2,295 
2,175 
2,607' 
2,949_^^ 
2.441 + 
2.488 



6l2 

+ 

-0 
+ 



2.99512 
3.04212 

3.14812 
2.80812 



1.871 



+ 
-0 

^t°o 



2.360 
2.627 
3.230 



+0 

+ 
-0 

=l°o 

^t°o 
«+o 



2.615 
2.444' 
2.948 
2.251 + 
3.182 
2.610 
2.773 
2.782 



2I2 
0I2 
3 + 



't2 
2.59712 
3.04212 

2.89812 
2.I75I2 
2.957I2 

n.5Q + 



3.059 
2.193 



2.898 
2.771 



028 
026 
082 
075 
050 
037 
030 
.029 
025 
.024 
045 
042 
072 
058 
034 
.032 
031 
028 
025 
024 
039 
036 
033 
.031 
014 
.014 
032 
.027 
060 
.056 
045 
043 
.013 
.013 
.046 
.039 
036 
.033 
008 
.008 
014 
.013 
031 
.030 
007 
.007 
021 
.019 
023 
023 
022 
021 
079 
060 
015 
.015 
010 
010 
009 
009 
024 
.024 
118 
.103 
078 
069 
028 
027 
.102 
.103 
.055 
.052 
.067 
.068 
.044 
.041 
.060 
.055 



1.064]^'^ 

I.O52I2 
, + 
-0 
I.227I2 
I.774I2 

1.60012 

0.690 ' " 



1.772^ 



2.800^ 



5±2 



oi2 
+0 
-0 
1+0 



0.931 

2.521_ 
,+0 
-0 
I.344I2 
O.921I2 

-^t2 
0.88412 



0.812^ 



0.485^ 



1.761 



+ 
-0 

7+0 
-0 



8l2 



.451 
.451 
.313 
.313 
347 
347 
.60S 
.608 
.484 
.484 
346 
346 
.229 
.229 
.559 
.559 
.297 
.297 
.292 
.292 
.318 
318 
511 
511 
385 
385 
522 
522 
981 
981 
474 
474 
.238 
.238 
.442 
.442 
991 
991 
186 
.186 
.369 
.369 
.575 
575 
225 
225 
944 
944 
214 
214 
319 
319 
449 
449 
127 
127 
230 
230 
479 
479 
466 
466 
.163 
.163 
.508 
.508 
.278 
.278 
.116 
.116 
.381 
.381 
.377 
.377 
753 
753 
.877 
.877 



9 


.3912 


.49 
.41 





.334+2; 


072 
.070 


2. 


.489 + 2 


.028 
.027 


1. 


74012' 


.425 
.425 


9 


.6612 


.99 
.85 





.59712; 


196 

,174 


3. 


.013^2 


.084 
.074 


0. 


681 + 2' 


.259 
.259 


6 


.6712 


.30 
.30 





,498 + 2' 


095 
.089 


2. 




.047 
.044 


1. 


018+2' 


.315 
.315 


8 


.-12 


.39 
.38 





,395 + 2' 


084 
.082 


2. 


,282 + 2' 


.030 
.029 


1. 


,95512 


.578 
.578 


10 


.3712 


.54 
.52 





,426+2' 


085 
.081 


2. 


.158+2 


.027 
.026 


1. 


.75312 


.430 
.430 


11 


' = 3lJ 


.24 
.94 





,25912: 


130 
.127 


2. 


,55012 


.048 
.044 


0. 


899+2 


.285 
.285 


6 


,6712 


36 
.35 





606 + 2' 


122 
.114 


2. 


,927+2 


.068 
.058 


0. 


,658 + 2 


216 
.216 


8 


57+0 


37 
.37 





411 + 0. 


077 
.075 


2. 


,40312 


.033 
.032 


1. 


,768 + 2 


513 
.513 


6 


,9312 


27 
.27 





369 + 2' 


075 
.072 


2. 


,44912 


.030 
.029 


0. 


,85712 


255 
.255 


7, 


.9612 


26 
.25 





,321 + 2' 


054 
.054 


2. 


93612 


.025 
.024 


0. 


923 + 2 


254 
.254 


10 


32 + 


.56 
.55 





,48112: 


099 
.095 


2. 


988 + 2 


.041 
.039 


0. 


,87112 


265 
.265 


7 


,9 + 0. 

■'■'-a 


32 
.31 





,47312; 


067 

.065 


3. 


,09512 


.032 
.031 


1. 


.51012 


.441 
.441 


6 


21+0 


09 
.09 





,45612; 


028 
.026 


2. 


,78412 


.014 
.014 


1. 


,169+2 


.367 
.367 


10 


39+0 


92 
.83 





,13712; 


116 
,111 


1. 


,85812 


.033 
.027 


1. 


,6isl2 


479 
.479 


6 


04+0 


.56 
.46 





,352 + 2' 


171 
,151 


2. 


.028 + 2' 


.061 
.056 


2. 


808 + 2 


989 
.989 


9 


,8512 


74 
.67 





,25812: 


103 
,102 


2. 


,625 + 2 


.044 
.042 


1. 


,370+2 


409 
.409 


7, 


,4612 


09 
.08 





,45812; 


025 
.025 


3. 


,234+2 


.007 
.007 


0. 


,605 + 2 


209 
.209 


13 


.10+ + 


.23 
.06 





.182+0. 


145 
.135 


2, 


.31712 


.050 
.041 


1, 


.15212 


.384 
.384 


5 


.8212 


.19 
.16 





.466+2' 


063 
.061 


2, 


.589+2 


.032 
.030 


2, 


.51212 


.938 
.938 


10 


8lt°0 


.12 
.12 





.28512: 


018 
.018 


3, 


.152+2 


.008 
.008 


0, 


.616+2 


.156 
.156 


13 


.3612 


.39 
.38 





.50412; 


044 

.043 


2 


.566+2 


.015 
.015 


0, 


.97212 


.299 
.299 


7, 


.0812 


.27 
.26 





.51612; 


078 
.076 


2. 


420+2 


.032 
.030 


1. 


.469+2 


513 
.513 


7, 


.60I2 


.07 
.07 





.45112; 


015 
.015 


2. 


,90412 


.007 
.007 


0. 


.84712 


205 
.205 


6 


.2812 


.15 
.15 





.54912; 


046 
,044 


2. 


233+2 


,024 
.024 


2. 


.25412 


.845 
.845 


10 


.5012 


.33 
.34 





.386+2' 


058 
.057 


3. 


112 + 


.024 
.024 


0. 


.678+2 


.180 
.180 


8 




.27 
.28 





.434+2: 


057 
.054 


2 


579 + 
■"''-0 


.023 
.021 


1. 


.18612 


282 
.282 


7, 


47+0 


.96 
.77 





.170+2' 


181 
.156 


2. 


750 + 


.080 
.059 


0. 


893 + 2 


437 
.437 


16 


72+0 
72_o 


.47 
.48 





.353+2' 


048 
.048 


2. 


699 + 2 


.017 
.017 


0. 


,372 + 2 


097 
.097 


11 


21+0 


.19 
.19 





.334I2: 


031 
.030 


2. 


,,543 + 2 


Oil 
.010 


0. 


755 + 


196 
.196 


6 


92 + 


.06 
.06 







.020 


2. 


986 + 2 


.009 


1. 


,59212 


.433 


6 


92 + 


.17 
.17 





,246 + 2' 


051 
.049 


2. 


R55 + 


.024 
.024 


1. 


,59612 


426 
.426 


6 


,ioi2 


.79 
.66 





,77312: 


275 
.239 


2. 


,176+2 


.119 
.105 


2. 


214+! 


135 
.135 


5 


,2612 


.28 
.27 





364+2' 


117 

.117 


2. 


,919+2 


.074 
.070 


0. 


,634+2 


.461 
.461 


7, 


.9012 


.26 
.25 





534+0. 


060 
.060 


3. 


00312 


.027 
.027 


0, 


83512 


.248 
.248 


6 


.231 + 


.26 
.92 





,190+2' 


217 
.180 


2. 


229I2 


.098 
.098 


2. 


,099+2 


.965 
.965 


9 


i«t°o 


.79 
.66 





.514+2: 


134 
,126 


2. 


622 + 2 


.056 
.052 


1. 


094+2' 


.353 
.353 


5 


.5512 


.23 
.22 





.663+2' 


115 
.106 


2. 


,878 + 2 


.064 
.066 


0. 


988 + 2 


.362 
.362 


8 


.81I2 


.64 
.58 





,418 + 2' 


106 
.104 


2. 


,71412 


.046 
.043 


2. 


,469 + 2 


682 
.682 


8 


.3312 


.67 
.61 





,46712: 


133 
.125 


2, 


,74212 


.060 
.055 


2, 


.50412 


.745 
.745 



tD 
CO 

CD 



Units are XIO 15 counts .s 
^ Using Njj = 12.1 X 10^0 
Using Njj = 18.7 X 10^0 



